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Dynamics of the singlet-triplet crystalline electric field (CEF) system at finite temperatures 
is discussed by use of the non-crossing approximation. Even though the Kondo temperature 
is smaller than excitation energy to the CEF triplet, the Kondo effect appears at tempera- 
tures higher than the CEF splitting, and accordingly only quasi-elastic peak is found in the 
magnetic spectra. On the other hand, at lower temperatures the CEF splitting suppresses the 
Kondo effect and inelastic peak develops. The broad quasi-elastic neutron scattering spectra 
observed in PrFe4Pi2 at temperatures higher than the quadrupole order correspond to the 
parameter range where the CEF splittings are unimportant. 

KEYWORDS: Kondo effect, dynamical susceptibility, inelastic neutron scattering, non-crossing 
approximation (NCA), filled skutterudite, PrFe4Pi2, PrOs4Sbi2 



1. Introduction 

Heavy fermion states have been observed in some Pr skutterudites such as PrOs4Sbi2 and 
PrFe4Pi2- PrOs4Sbi2 is the first heavy fermion superconductor discovered in Pr compounds 1 
and, in addition, has quadrupole ordered phase under high magnetic field. 2 PrFe4Pi2 exhibits 
antiferro-quadrupole order at low temperatures. 3 In the inelastic neutron scattering experi- 
ment of PrFe4Pi2, any CEF excitation peak cannot be detected at temperatures higher than 
the phase transition, while inelastic peak develops at lower temperatures. 4 On the other hand, 
the CEF excitations are clearly seen in PrOs4Sbi2- 5 ~ 7 Taking account of the fact that only 
PrFe4Pi2 exhibits the Kondo effect in Pr skutterudite series, 8 we consider that the heavy 
mass in PrFe4Pi2 is caused by the Kondo effect, while another mechanism may be relevant to 
PrOs 4 Sbi 2 . 

CEF level structures are essential to the properties of 4/ electrons. The CEF singlet ground 
state and low lying excited triplet states make clear the origin of the high field ordered phase 
in PrOs4Sbi2- 9 ' 10 The ordered phase in PrFe4Pi2 is also ascribed to the quadrupole moment 
of the CEF singlet-triplet system. 11 Difference of the actual wave functions of the CEF triplet 
states should create the diversity of physical properties in Pr skutterudite. 

According to the band calculation, main conduction band in Pr skutterudites is formed 
by the 2p orbitals of pnictogens with a u symmetry, which does not have the orbital degrees 
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of freedom. 12 Considering exchange interaction between 4/ 2 states and a u orbital, we have 
revealed the condition for the occurrence of the Kondo effect. 13 Then the fact that the Kondo 
effect appears only in PrFe4Pi2 is reasonably understood by the difference between the triplet 
wave functions. 

When the Kondo effect occurs in the singlet-triplet system, the Kondo effect competes 
with the CEF singlet. It is expected that the competition affects the physical quantities, 
especially their temperature dependence. The Kondo effect in the CEF singlet ground state 
has been studied theoretically by several authors in relation to UPai2Si2. 14 ~ 16 However, its 
dynamics such as dynamical susceptibility has not been derived in spite of its significance. 

In this paper, we derive dynamics of the CEF singlet-triplet system, and clarify tempera- 
ture dependence of dynamical quantities. The competition between the CEF singlet and the 
Kondo effect is the main issue. Dynamics is derived by the non-crossing approximation (NCA) 
in §2. Numerical results by the NCA and by the numerical renormalization group are given 
in §3. Then we discuss application to the real materials. We summarize the results in the last 
section. Details of the NCA equations for an actual computation are given in Appendices A 
and B. 

2. Application of the NCA to the Singlet-Triplet System 

The NCA has been developed as a solver of the Anderson model. 17 ' 18 The NCA enables 
us to compute dynamical quantities at finite temperature as well as static ones. An equivalent 
approximation has been applied to the Coqblin-Schrieffer model with CEF splittings, and 
temperature dependence of dynamical susceptibility has been discussed. 19 We shall newly 
apply the NCA to the localized 4/ 2 system with exchange interactions. 

2.1 Singlet-triplet Kondo model 

We have derived the effective exchange interaction in the singlet-triplet system coupled 
with the a u conduction spin s c . 13 The interaction is written in terms of operator X t and X s , 
or pseudo-spins Si = (X t + X s ) /2 and S 2 = (X t - X s ) /2 as follows: 

H s . t = e t P t + e s P s + (7 t X* + I S X S ) ■ s c 

= A CEF S 1 -S 2 + (J 1 S 1 + J 2 S 2 )-s c , (1) 

where P s = — Si ■ S 2 + 1/4 and P t = Si ■ S 2 + 3/4 are the projection operator onto singlet 
and triplet states, respectively. Ji :2 = I t ± I s and Aqef = £t — £s is the CEF splitting. 
We have chosen the origin of energy so that 3e t + e s = 0. In the second order perturbation 
on hybridization, the coupling constant Ji becomes ferromagnetic and J 2 antiferromagnetic, 
provided the triplet wave functions mainly consist of ^4) in the point group Oh- On the other 
hand, both exchange interactions can be negligible if the triplet is almost composed of ^5) in 
O h . 
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Fig. 1. Operators acting on the singlet-triplet system. 




Fig. 2. Diagrammatical representation of the equation for the resolvent and self-energy in the NCA. 
Solid, dashed and wavy line denote conduction electron, resolvent and interaction, respectively. 
Single and double line represent bare and renormalized one, respectively. 



The operator X act on the triplet states, and X s connect the singlet and the triplet as 
shown in Fig. 1. In applying the NCA to the singlet-triplet Kondo model, it is convenient to 
work in the form of X t and X s because singlet and triplet are properly distinguished from 
each other. 

2.2 Integral equations for resolvents and effective interactions 

We derive dynamics of the singlet-triplet system adopting the NCA. In applying the NCA 
to the exchange interactions, we introduce a fictitious 4/ 1 intermediate state with negligible 
population. 18 We could alternatively introduce a fictitious 4/ 3 interdediate state instead of 
the 4/ 1 state. 

We introduce resolvent R s (z) for the singlet state and Rt(z) for the triplet. The effects of 
the interactions are taken account of as the self-energy T, a (z). Each resolvent R a (z) is given 
by 

R a (z) = [ z -e a -Z a (z)]-\ (2) 

where a denotes configurations of 4/ 2 states. The NCA determines Ti a (z) in a self-consistent 
fashion. Figure 2 shows diagrammatical representation of the equation for the resolvent and 
self-energy in the NCA. 

We obtain the renormalized exchanges It(z) and I s (z), which are modified from the bare 
ones It and I s . In addition, an effective potential Kt(z) for the triplet and K s (z) for the singlet 
are generated by higher-order exchange scatterings. The self-energy is given in terms of the 



3/16 



Full Paper 




Fig. 3. Diagrammatical representations of the equations for the effective interactions in the NCA. 



effective potentials by the NCA integral equation: 

£ Q (z) = -2 J dep c {e)[l - f(e)]K a (z - e). (3) 

In order to derive equations for the renormalized interactions I a (z) and K a (z), we divide the 
operator X s into two parts: 



X s = P t X s P s + P s X s P t . 



(4) 



The first term operates on the singlet state, and second one to the triplet states. Corre- 
spondingly, we define effective interactions Is ts \z) and 1$ (z) for each part. The effective 
interactions are determined by equations illustrated in Fig. 3. The simultaneous equations for 
the effective interactions are given in the matrix form by 
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(5) 



dep c (e)/(e)lk(z + e), (6) 

which appears as a bubble diagram in Fig. 3. In deriving eq. (5), we have used the following 
properties: 

) = (2P t <r u + X 1 ■ a) 

act' ? 

) = (P t <r° + 3P s a° + X t -a) 

ai<Ti 

^ (-^aai ' °"°-io-')(^a ia ' ' °Wi) = 2Pt-X'^ a ,P s • <T CT(T / , 

Q1CT1 
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^2( X a ai ■ T<Ti*')( X a 1 a' ' a <™l ) = 2P s X aa' P t " ° oa' , (7) 
aitri 

where <r° denotes a unit operator for the conduction spin. Solving eq. (5) for the renormalized 
interactions, we obtain 

~ = 4(2J t - / s 2 n s ) ci 

* (2-/ t n t )(4 + 4/ t n t -3/|n t n s ) a &' 

J = f(ts) = f(st) = = ^2 

s ~ s s (4 + 4/ t n t - 3/|n t n s ) ~ b ' 

~ 4i t 2 n t + 2/ 2 n s - 3/ t / 2 n t n s _cs 

* (2 - j t n t )(4 + 4/ t n t - 3/ s 2 n t n s ) ~ a &' 

^ = 3/ 2 n t _C4 

4 + 4/ t n t - 3/ 2 n t n s b ' i; 

where the simplifying notation a = 2 — / t n t , 6 = 4 + 4/ t IT t — 3/ 2 n t n s , etc. will be utilized in 
Appendix A. The effective interactions I a (z) are required only by the vertex corrections for 
the dynamical magnetic susceptibility. 

2.3 T -matrix 

In the NCA for the exchange interactions, the impurity T-matrix is computed as follows: 

T(ie n ) = -±- [ ^e-P z J2Ka(z)R a (z + ie n ), (9) 
z f Jc 2m V 

where e n = (2n + 1)ttT is the Matsubara frequency of fermions and the contour C encircles 
all singularities of the integrand counter-clockwise. We utilize the spectral intensities 

r/ a (uj) = --ImR a (uj + iS), (10) 

7T 

where <5 is positive infinitesimal. We further define rjcPfa) and £^\uj) for K a (z), and r]a {u) 
and £a^(o>) f° r Ia{ z ) m a similar fashion. The spectral function ^ a {uj) should be actually 
computed by another set of equations to avoid difficulty of the Boltzmann factor at low 
temperatures. We shall describe the details in Appendix A. Performing analytic continuation 
ie n — > uj + id to real frequencies in eq. (9), we obtain 

— ImT( W + tf) = £ / de[^ K He) Va (e + u) + V i K \e)Ue + ^}- (12) 
Electrical conductivity is derived from lmT(oo) as described later. 

2.4 Dynamical magnetic susceptibility 

The dynamical magnetic susceptibility is given by 

X{iv m ) = {gj^B) 2 ^2 (a'\J z \f3)(f3'\J z \a)x(aa',f3f3';iu m ), (13) 

ota'PP' 
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Fig. 4. Diagrammatical representation of the dynamical magnetic susceptibility. 



Table I. Combinations (pta' , (3/3') and relevant quantities for the dynamical magnetic susceptibility. 



act', (3(3' 


(a'\J z \P){(3'\J z \a) 


A(aa', (3(3'; z,z + iv m ) 


t+t+, t+t+ 


|(t + |J,|t+)| 2 


A t (z,z + iv m ) 


t-t-, t-t- 


l(t-|J,|t-)| 2 


A t (z,z + iu m ) 


s s, tO tO 


|(tO|J z js)| 2 


A sp (z,z + iv m ) 


tO tO, s s 


|(tO|J z |s)| 2 


A sp ,(z, z + iu m ) = A sp (z + iu m , z) 


s tO, s tO 


<tO|J z |s) 2 


A sc (z, z + iv m ) 


tO s, tO s 


<s|J,|tO} 2 


A sc /(z, z + iu m ) = A sc (z + iu m , z) 



X(aa,(3(3';iv m ) = -^- f ^-.e l3z A(aa ', (3(3'; z, z + iv m ), 
Zf Jc 2m 



(14) 



where v m = 2mirT is the boson Matsubara frequency, and A(aa' , (3(3'; z,z + iv m ) is the vertex 
function to be explained below. The operator J z has non-zero matrix elements for 

(t ± \J z \t±) = ±M t , (s|J z |tO) = M s . (15) 

Although M t and M s depend on the triplet wave functions in the point group T^, 10 we take 
these matrix elements as unity in the following numerical calculations. It is straightforward 
to reinstate the actual values for particular wave functions. 

The vertex part A(aa' , (3(3'; z, z+iv m ) is diagrammatically represented in Fig. 4. It contains 
infinite pairs of conduction electrons and holes through the Bethe-Salpeter-type equations in 
the NCA. There are six combinations (act', (3(3') that contribute to the dynamical susceptibil- 
ity. A combination is symbolically represented by A hereafter. Only three kinds of vertices are 
independent (see Table I). We note that (t + | J z \t+) and (t — |J z |t— ) are completely decoupled 
owing to cancellation between the processes through the state |t0). Figure 5 shows diagram- 
matical representation of the equations for the vertices. The equations are analytically given 
by 

A t (z,z') = R t (z)R t (z') + R t (z)R t (z')2 f dep c (e)[l - /(e)] J de' p c (e')f(e') 

A t (z-e + e',z' -e + e'), (16) 

A sp (z,z') = R s (z)R t (z') + R s (z)R t (z')2 j dep c (e)[l - /(e)] J de' p c (e')f(e') 



-h{z - e)h(z' - e) + K t (z - e)K t (z' - e) 
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Fig. 5. Diagrammatical representation of the equations for the dynamical magnetic susceptibility. 



x [K s (z - e)K t {z' - e)A sp (z -e + e',z'-e + e') 

+~I a (z - e)h(z' - e)A sc (z - e + e',z'-e + e% (17) 

A sc (z,z') = R t (z)R s (z')2 j dep c (e)[l -/(e)] J de' p c (e')f(e') 
x [K t (z - e)k s (z' - e)A sc (z - e + e',z r - e + e') 

+ -J s {z - e)Uz' - e)A sp (z -e + e',z'-e + e% (18) 

The other vertices A 8p >(z, z') and A sc '(z,z') satisfy the equations with replacement Rt <-> Rs, 
K t <-> K s , A sp — * A sp ', and A sc —* A sc / in eqs. (17) and (18), respectively. 

After analytic continuation %v m — > uj + id to the real axis in eq. (14), we obtain 

Im XA (u, + iS) = (1 - e-^)— / —e-* 

Zf J_ 00 2ir 

x Re[A A (e - iS, e + u + iS) - A A (e + id, e + to + i<S)]. (19) 

The static susceptibility xa(0) is given by 

1 f°° df 

Xa(0) = — / -e-^ImA A (e + iff, 6 + iff). (20) 

In numerical calculations, it is difficult to perform the integration including the Boltzmann 
factor in eqs. (19) and (20). We avoid the difficulty through transforming eqs. (16), (17) and 
(18) into equivalent and convenient forms (see Appendix B). 

3. Numerical Results for Resistivity and Susceptibility 

The singlet-triplet Kondo model has two characteristic energy scales, i.e., the Kondo tem- 
perature Tk and renormalized CEF splitting Acef- Provided the CEF singlet has lower energy 
than the CEF triplet, the Kondo effect and the CEF singlet compete with each other. We 
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Fig. 6. J\-Ji phase diagram of the ground state with Acef = 0. 



shall examine how the competition affects the temperature dependence of dynamical quanti- 
ties. We first show the phase diagram of the ground state by the numerical renormalization 
group (NRG), and next dynamical quantities obtained by the NCA. 

3.1 Numerical renormalization group analysis 

The singlet-triplet Kondo model has four fixed points, i.e., CEF singlet, CEF triplet, 
doublet and quartet. They have residual entropies In 1 = 0, In 3, In 2 and In 4, respectively. 
Some explicit results for temperature dependence of the entropy by the NRG have been 
presented in ref. 13. Figure 6 shows the J1-J2 phase diagram of the ground state of the 
singlet-triplet Kondo model with the bare CEF splitting Acef = 0. Note the symmetry 
against interchange of J\ and J2 . The effective CEF splitting becomes finite by renormalization 
effect. In the second order perturbation theory, self-energies of the triplet and singlet states 
are written as 

*?\z)~- l -D P l{2lt+I^ 

^(z)~-\ Dp lll (21) 

(2) 

at T = and \z\ <C D. Then we obtain the renormalized CEF splitting A C ^ F in the second 
order perturbation theory as 

4ef = -DplVs ~ ft = -D P 2 c JiJ 2 . (22) 

It is now clear that with J1J2 < 0, the second order exchange stabilize the CEF singlet. 
On the other hand, higher order exchanges cause the Kondo effect, which makes the CEF 
splitting obscure. The competition between the Kondo effect and the CEF effect depends on 
their characteristic energy scales. We find that renormalization of the CEF splitting is larger 
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than the Kondo temperature in the calculated parameter range. 

In computing dynamical quantities, we take J 2 = to minimize renormalization of the 
CEF splitting. Even if J2 equal zero, pseudo-spin S2 would interact with conduction elec- 
trons indirectly through another pseudo-spin Si. The singlet-triplet Kondo model with the 
additional condition Acef = is reduced to the Kondo model where the pseudo-spin S2 is 
decoupled. 

3.2 Electrical resistivity 

We discuss how the CEF singlet influences the Kondo effect, comparing with the behavior 
in the CEF triplet ground state. We take a rectangular model with the band width 2D for 
conduction electrons 



where A is a constant and /(e) is the Fermi distribution function. Figure 7 shows temperature 
dependence of the electrical resistivity p{T) = l/a(T) with J\p c = 0.2. We take Acef rang- 
ing from — 40K (triplet) to 40K (singlet). In the case of the triplet ground state, resistivity 
monotonously increases as temperature decreases, corresponding to development of the Kondo 
resonance peak. The CEF splitting only causes a slight reduction of the slope in comparison 
with the quartet case Acef = 0. On the other hand, the CEF singlet leads to suppression of 
the enhancement at temperatures below about the third of Acef- The suppression is due to 
a pseudo-gap of about 2Acef in ImT(u;) around Fermi level. Note that the peak temperature 
of p(T) is substantially smaller than Acef- 



p c {e) = 9(D - \e\)/2D 



(23) 



where 6(e) is a step function. We tentatively assume D = 10 4 K and take Kelvin as the unit 
of energy. 

We obtain the electrical conductibity cr(T) by 




(24) 
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Fig. 7. The electrical resistivity as a function of temperature for several values of Acef- 
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3.3 Dynamical magnetic susceptibility 

We show typical behaviors of magnetic spectra Imx(w) /(l — e~^) both in the CEF singlet 
and triplet ground states. Figure 8 shows the temperature dependence of magnetic spectra. 
In the case of the triplet ground state, quasi-elastic peak develops sharply with decreasing 
temperature. Excitation to the CEF singlet is also found at to ~ 30K. In the case of the CEF 
singlet ground state, quasi-elastic peak is seen at temperatures higher than Acef> an d an 
inelastic peak develops at lower temperatures. The temperature dependence is ascribed to the 
competition between the CEF singlet and the Kondo effect. Namely, the behavior shows a 
change from the dominance of the Kondo state at higher temperatures to the CEF singlet at 
lower temperatures. 




Fig. 8. Temperature dependence of the magnetic spectral intensities. J\p c = 0.2, J 2 p c = 0, and CEF 
splittings Acef are — 40K (triplet) in (a) and 20K (singlet) in (b). 



The magnetic relaxation rate T measured by NMR is given by 

lim ^1 = 242) . (25) 

If Imx(u;) can be approximated by a Lorentzian, T corresponds to the half-width of a quasi- 
elastic peak. Figure 9 shows the temperature dependence of T for several values of Acef- In 
the temperature region where the Kondo effect occurs, F behaves as T a with the exponent a 
being smaller than unity. In the case Acef > 0, the crossover to the CEF singlet from the 
Kondo effect causes an increase of V as temperature decreases. It results from the following 
reason: Imx(w)/w |^_^o is reduced at lower temperatures while x(0) is constant due to the van 
Vleck susceptibility in eq. (25). The decrease in the case of Acef = 20K and 40K below about 
5K is caused by the inaccuracy of the NCA at low frequencies, which gives fictitious increase 
of Imx(w)/u;. 21 It is known that half-width of Imx(w)/w at zero temperature corresponds to 
the Kondo temperature Tk. We estimate as Tk ~ IK for Acef = from Fig. 9. 

Next we examine the Kondo temperature Tk dependence. Figures 8(b) and 10 show mag- 
netic spectra with different values for coupling constant Ji but with renormalized CEF split- 
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Fig. 9. The magnetic relaxation rate T as a function of temperature for several values of Acef- 



ting Acef fixed. These parameters satisfy the condition Tk < Acef, i-e., no residual entropy. 
In a small coupling case J\p c = 0.18, inelastic peak is seen more clearly as compared with 
the case of J\p c = 0.20. On the other hand, a large coupling lowers the crossover temperature 
from the Kondo effect to the CEF singlet. Inelastic peak does not arise at all temperatures 
calculated in Fig. 10(b). The apparent inelastic peak at u ~ 6K is actually a quasi-elastic 
scattering whose maximum is shifted from zero by the Bose factor. The crossover temperature 
depends both on Tk and Acef- 




Fig. 10. Temperature dependence of the magnetic spectral intensities with different values for cou- 
pling constant J\ but with almost the same value for the renormalized CEF splitting Acef ~ 10K. 
Jip c = 0.18, Acef = 16K and (b) J lPc = 0.22, A C ef = 22K. 



3.4 Application to Pr skutterudites 

We discuss the origin of diversity of the magnetic spectra in PrOs4Sbi2 and PrFe4Pi2 
on the basis of results obtained for the singlet-triplet Kondo model. In the case of the CEF 
level structure for PrOs4Sbi2, the 4/ 2 configuration has negligible antiferromagnetic exchange 
interaction with conduction spins because the low lying triplet consists mainly of in Oh- 13 
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We accordingly expect clear CEF excitation, which has actually been observed. 

In the inelastic neutron scattering experiment of PrFe4Pi2, on the other hand, only broad 
quasi-elastic peak is found at temperatures higher than the phase transition. We ascribe the 
quasi-elastic scattering to the Kondo effect, the occurrence of which is consistent with the 
temperature dependence of the electrical resistivity. 8 Inelastic peaks are detected only at 
lower temperatures. The observed inelastic peaks involve a suppression of the Kondo effect. 
The suppression may be due to either the quadrupole order or formation of the CEF singlet. 
We conjecture that the former origin is likely for PrFe4Pi2, in view of recently proposed CEF 
level scheme, i.e., nearly degenerate singlet-triplet levels. 20 The CEF splitting in that scheme 
is not large as compared with the Kondo temperature. Thus the inelastic neutron scattering 
experiment can resolve the levels only when the quadrupole order breaks the Kondo effect. 
Alternatively the level splitting of the triplet by internal fields is observed. 

4. Summary 



We have derived dynamics of the CEF singlet-triplet system coupled with conduction 
spins. Explicit calculation is performed with use of the NCA. The physical quantities are 
characterized by two energy scales, i.e., the Kondo temperature 7k and renormalized CEF 
splitting Acef, which is modified from bare one Acef due to exchange interactions. 

Provided the condition Ik < Acef is satisfied, ground state is the CEF singlet. Although 
4/ 2 configuration forms the CEF singlet at low temperatures, the Kondo coupling can be 
present at temperatures higher than Acef- Correspondingly, the electrical resistivity is en- 
hanced with decreasing temperature, while it is suppressed by the CEF splitting eventually. 
The competition between the Kondo effect and the CEF singlet also leads to peculiar depen- 
dence of the magnetic spectra on temperature. As shown in Fig. 8(b), the inelastic peak cannot 
be resolved from the broad quasi-elastic peak at T > Acef- The inelastic peak, however, may 
be identified as temperature decreases to T < Acef- 

We have discussed the difference between PrFe4Pi2 and PrOs4Sbi2 in terms of different 
wave functions of the CEF triplet states. This explains the contrasting behavior observed in the 
inelastic neutron scattering experiment. In PrFe4Pi2, the first excited CEF triplet is composed 
mainly of T4 in Oh, and antiferromagnetic coupling between the CEF states and conduction 
spins becomes significant. Accordingly, PrFe4Pi2 can exhibit the Kondo effect even with the 
CEF singlet ground state, and the CEF excitation is suppressed. Since our calculation does not 
include intersite interactions, however, we cannot exclude the possibility that the quadrupole 
order is responsible for suppression of the Kondo effect and the development of the inelastic 
feature at low temperatures. 
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Appendix A: Equations for the spectral functions 

Physical quantities are conveniently represented by two types of spectral functions, i.e., 
rj a (io) and £ Q (u;). Although £ a (cj) is related to 7] a (uj) analytically, it is difficult to com- 
pute directly due to the Boltzmann factor at low temperatures. Following ref. 21, we ob- 
tain £ a (a;) by numerical iteration as follows. We introduce an operator by R(uj) = 
— Zj l e~ l3u) lmR(uj + i5)/ir. Operating on eqs. (3) and (6), we obtain 

a a (oo) = P«) Ea ( w ) = -2 J dep c (e)f(e)& K \oo - e), (Al) 

n a (u) = P®n a (u) = J de Pc (e)[l - /(e)]£ a (c; + e). (A-2) 
Here £, a (oj) are related with the corresponding resolvent as 

CaH = V^R a {u) = \R a (u + iS)\ 2 a a (u). (A-3) 

Analogously, we obtain 

aU) = p(0/ t = J^i![(p(0 a *)6* Cl + a {P®b*)a + ab(V^ Cl )}, 
c i 



C2 



e W _ p(0 ^ = _^tP [(73 «) a * )6 * C3 + a (p(0 6 *) C3 + a 6(P (f) c 3 )], 
l c 3r 



= 7?«)^ s = _L^[ (77 (5) 6 *) C4 + b(V^c 4 )}, (A-4) 

l c 4| 

where we have used the notations in eq. (8) and 

V®a* = / t 7r t , 

p(Oft* = _47 t7 r t + 3/ s 2 (vr t n s + n t ir B ), 

V^c 3 = 4/ t V t + 2/ s 2 vr s - 3/ t / s 2 (^ t n s + n t vr s ), 

V®c A = 3/ s V t . (A-5) 
The norms are determined by the following sum rule: 

/oo 
dw&M = !■ (A-6) 

Normalization by eq. (A-6) is equivalent to determination of the partition function Zf of / 
electrons. 
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Appendix B: Equations for the magnetic susceptibility 

The formulae of Imx(to) in eq. (19) and x(0) in eq. (20) include integrations with the 
Boltzmann factor. To avoid this factor, we rewrite the equations following ref. 22. For each 
frequency u we introduce the following simplifying notations: 

Q\(e) = A A (e + iS, e + u + iS), 

Q\(e) = A A (e - i5, e + u + i5), 

R aa '(e) = R a {e + i5)R a >(e + u + i5), 

R aa >(e) = R a (e-i5)R a ,(e + u> + i5). (B-l) 

Iaa>(c), laa'i^)-, K aa i(e) and K aa >(e) are defined with use of I a {z) or K a (z) in a similar 
manner. We introduce an operator V by VQ{e) = Zj 1 e' l3e [Q(e) — Q(e)]. In terms of V, eq. 
(19) is represented by 

/°° eh 
—ReVQx(e), (B-2) 

without the Boltzmann factor. Operating V on eqs. (16), (17) and (18), we obtain 
VQ t (e)=PR tt (e) 

+ [VR tt (e)]2 J de'p c (e')[l-f(e')} J de"p c (e")f(e")[±I tt (e-e') + K tt (e-e')]Q t (e-e' + e") 

+ R tt (e)2 J de' Pc (e')f(e') J de" p c (e")f(e")[±VI tt (e - e') + VK tt (e - e')]Q t (e - e' + e") 

+ R tt (e)2 J de' Pc (e')f(e') J de"p c (e")[l - f(e")][\l tt (e - e') + K tt (e - e')][VQ t (e - e' + e")], 

(B-3) 

VQ sp {e) =VR st (e) 

+ [PR Bt (e)]2 J de>p c (e')[l-f(e')} J de" p c (e")f(e")K st (e - e')Q sp (e - e' + e") 

+ R st (e)2 J de>p c (e')f(e>) J de" p c (e")f(e")[VK st (e - e')]Q sp (e - e> + e") 

+ i? st (e)2 J de'p c (e')f(e') J de"p c (e")[l - f(e")]K st (e - e')[VQ sp (e - e' + e")] 

+ [VR st (e)]2 J de'p c (e')[l-f(e')} J de" p c (e")f(e")I ss (e - e')Q sc (e - e' + e") 

+ R st (e)2 J de' Pc (e')f(e') J de"p c (e")f(e")[VI ss (e-e')]Q sc (e-e' + e") 

+ £ st (e)2 J de' Pc (e')f(e') J de"p c (e")[l - f(e")]I ss (e - e')[VQ sc (e - e' + e'% (B-4) 

VQ sc (e) 
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= [VR ts (e)}2 J de'p c (e')[l-f(e')} J de" p c (e")f(e")K ts (e - e')Q sc (e - e' + e") 

+ R ts (e)2 J de'p c (e')f(e') J de"p c (e")f(e")[PK ts (e-e')]Q sc (e-e' + e") 

+ R ts (e)2 J de' Pc (e')f(e') J de"p c (e")[l - f(e")]K ts (e - e')[VQ sc (e - e' + e")] 

+ [PR ts (e)]2 J de' Pc (e')[l-f(e')} J de" p c (e")f(e")I ss (e - e')Q sp (e - e> + e") 

+ R ts (e)2 J de' Pc (e')f(e') J de"p c (e")f(e")[VI ss (e-e')]Q sp (e-e' + e") 

+ R ts (e)2 J de'p c (e')f(e') J de" p c (e")[l - f(e")]I ss (e - e')[VQ sp (e - e' + e")]- (B-5) 

VR aa >{e) is given by known quantities as follows: 

VR aa ,(e) = 2iri£ a (e)R a/ {e + u + i6). (B-6) 

VI aa '{e) and VK aa /(e) are obtained by analogous equations. After Q 7 are obtained using eqs. 
(16), (17) and (18), VQ 1 are computed by numerical iterations of eqs. (B-3), (B-4) and (B-5). 

For the static susceptibility, we introduce another operator P by PQ(e) = Zj 1 e^ f3t ImQ(e) 
instead of V. The static susceptibility is given in terms of V as follows: 

r°° de - 

X 7 (0) = / -PQ 7 (e). (B-7) 

The equations for PQ 1 are obtained by replacing P, R aa , , I aa > , K aa , by P, R* aa , , I* a , , K* a , in 
eqs. (B-3), (B-4) and (B-5). We also obtain PR aa ,(e) by 

PRaa'(e) = -Tr[U^Ra'(e + i6) + R* a (e + iS)^(e)}. (B-8) 

The formulae for PI aa >(e) and PK aa >(e) are obtained by analogy of eq. (B-8). PQ 1 are 
numerically computed in a way similar to PQ-y. 
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